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Qh' Abstract. The observed distribution of globular cluster binary radio pulsars 

^ I in the eccentricity versus orbital period plane can be explained as a result of 

binary-single star interactions. Our numerical and analytical study hints that 
\^ ' the highest eccentricity binaries in clusters are likely to be from exchange and/or 

merger of a single star with a binary component, while the intermediate eccen- 
tricity systems are probably results of fly-by interactions. 

Q^. 1. Introduction 

2 ■ Observable parameters of neutron star (NS) binaries in globular clusters (GC), 

' such as spin period (Pspin), orbital period (Porb), eccentricity (e), projected 

. radial position in the cluster, companion mass etc. provide a valuable test-bed to 

examine the theoretical scenarios of formation and evolution. These parameters 
can even be used as a tracer of the past history of dynamical interactions of 
I the binary NSs in individual GCs. As orbital parameters of LMXBs are poorly 

. known, we concentrate entirely on binary radio pulsars whose orbital parameters 

can be determined accurately by pulsar timing. We study the distribution of 



o 



^f) • GC pulsars in e — Porfe plane as a result of interactions with single stars. 

o . 

G^ , 2. Observed data 

o 

So far 140 pulsars in 26 GC^ have been discovered, 74 are binaries, 59 are 
isolated and 7 without published timing solutions. PSR J2140— 2310B with 
poorly known orbital parameters is excluded from the present study. Majority 
' of the rest 73 binary pulsars are millisecond pulsars whose e Vs Port distri- 

bution is shown in Fig [TJ They can be categorized into three groups: I) 21 
pulsars with large eccentricity (1 > e > 0.01); II) 20 pulsars with moderate 
eccentricity (0.01 > e > 2 x 10~^); III) 32 pulsars with small eccentricity (e = 
in the database, but e = 3 x 10~ ^ in present w ork). The smallest measur 
able value of e follows the relation ( Phinnev 19921 ): Cmin = {6t)/{asin i/c) 

4TT^c6t/ [sin i (G(mp + mc)^/^) 



We take lUp (pulsar mass) = 1.4 Mq, 

rric (companion mass) = 0.35 Mq, i (inclination angle) = 60°, 6t (timing accuracy) 
1 usee to draw this "limit of timing sensitivity" line (the lower left corner 



^http://www.naic.edu/'^pfreire/GCpsr.html 
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of Fig. [T]). Observational selection effect may be influencing the distribution in 
the low Porb region (Camilo & Rasio 2005). 



3. Stellar interactions in globular clusters 

Due to high stellar densities, binary - single star interaction is significant in GCs 
resulting any one of among "fly-by" , "exchange" , "merger" or "ionization". 
Fly-by, exchange and merger interaction can produce eccentric binary millisec- 
ond pulsars when the interaction timescale is less than the binary age wheras 
formation scenario of millisecond pulsars suggests that their eccentricity should 
be very small, ~ 10"^ - 10"^ (Phinney 1992). We take the maximum age of the 
binaries in a GC to be equal to GC ages ('^ 10^'' yrs). Rasio &: Heggie (1995) 
obtained the expressions for fly- by time-scales for an initially circular binary as: 

tfiy = 4x 10"n4 it>ioPj/^e2/5 for e < 0.01 

tfiy = 2x IQi^n^ i7;ioPj/^ [- ln(e/4)]^'/3 for e > 0.01 

where is the number density (n) of single stars in units of lO'* pc~^ and vio is 
the velocity dispersion (v) in units of 10 km/sec in GCs; Port is the orbital period 
in days giving iy^j^ in years . The value of y i n / (in GC center) varies from 0.0024 
- 2.167 for different GCs (IWebbink 19851 ^. We grouped the GCs into SIX groups 
according to the values of vio/n^ and calculated tfiy with the mean values of 
^^10/^^4 for each group. In Fig. [H we plot the isochrones of fly-by encounters (tfiy) 
in the e — Porfe plane for all groups. If tjiy > 10^^ years for a particular binary, 
then it would not have interacted and it would preserve its original eccentricity. 
If tfiy < lO^'^ years for a particular binary, then it could be eccentric due to 
fly-by interactions which is the case for most of the eccentric GC binaries (Fig. 
[1]). Therefore, many GC pulsars' orbital eccentricities can be explained by fly- 
bys. However the local stellar densities at the pulsar positions (especially with 
positional offset from cluster cores) may result higher values of ^'lo/?^4 from the 
central values used here. Thus these pulsars might be in tfiy > 10^^ years region. 
Moreover, for pulsars with e > 0.1 (a majority of the group I pulsars), very close 
fly-bys are necessary if the initial binaries were circular (Camilo &: Rasio 2005). 
In these cases, exchange and merger events are relevant. But the three pulsars 
with 0.01 < e < 0.1 and 60 < Port < 256 d are p ossibly white dwarf c ores of 
red giant companions that overflowed Roche lobe ( Webbink et al. 19831 ). Such 



binaries would normally have the "relic" eccentricities ~ 10~^ (Phinney 1992). 
The above binary pulsars with their presently mildly high eccentricities, have 
undergone fly-by encounters with single stars, rather than exchange reactions, 
which would produce very high eccentricities e > 0.1. 

We used STARLAEg codes to perform numerical simulations of exchange 
and merger with different stellar parameters (GC parameters same as Terzan 
5, since it has the lowest value of the parameter vio/n^ where all classes of 
encounters can take place). In Fig. [21 the left y axis gives the final eccentricities 
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Figure 1. GC binary pulsars in the e — Porb plane with contours of tfiy — 
-j^gio yj.g (golid lines) and tfiy = 10^ yrs (dashed lines) for different values 
of viQ/n4 and contours of tgr — 10^° yrs (solid lines) and tgr = 10^ yrs 
(dashed lines) for binaries with mp — 1.4 Mq and rric — 0.35 Mq or 
0.16 Mq. The color scheme for tfiy and tgr contours are shown on the right. 
Individual pulsars are marked with same colors as vio /n^ values of their host 
GCs. Pulsars with projected positions inside the cluster core are marked 
with +, those outside the cluster core with x and the pulsars with unknown 
positions with o. Six binaries in Terzan 5 with e > 0.1 are marked with their 
names. 

while the right y axis gives the time scales of interactions. Moreover, we plot 
Porb,in of the initial binary along the top x-axis and Porbjin along the bottom 
X-axis. Porbjin is obtained from Porb.in putting A = in the relation ajj„ = 
[ainmanih / 17111712 (1 — A)] where mi and m2 are masses of the members of the 
initial binary, is the mass of the incoming star, rua and are masses of the 
members of the final binary, A is the fractional change of binary binding energy. 
For exchange, rua = rrii, rrih = rns; for merger, = mi, rub = 1712 + 7713. 
Six high eccentricity (e > 0.1) binaries in Terzan 5 are also shown (red in 
colors, symbols same as used in Fig [T|) here. It is clear from the scatter plots 
(Fig. [2|) that the final binaries will most probably have e > 0.1 if they undergo 
either exchange or merger events. If a pulsar lies in a region where tint > 10^'^ 
years in a particular plot, then that pulsar can not experience that interaction. 
Moreover, it the companion masses are much different from expected values 
(ms for exchange and m2 + for merger) then that interaction should also 
be excluded. Thus we conclude that PSR I (in Ter5) can be either result of 
exchange if m2 was either 0.16 Mq or 0.40 Mq or merger if 771-2 was 0.16 Mq, 
PSR J and X can be results of exchange if m2 was either 0.16 Mq or 0.40 Mq, 
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PSR Q, U and Z can be results of merger only if m2 was 0.40 Mq where in all 
cases we consider ni^ = 0.33 Mq (average mass of single stars in GCs). 
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Figure 2. Time scales (denoted by black '+') and final eccentricity distri- 
butions (colored scatter plots) with initial and final orbital periods for ex- 
change (purple points on the left panel) and merger (green points on the right 
panel) interactions with different stellar parameters. Vertical orange lines give 
boundaries of orbital periods where interaction time scales < 10^*^ yrs. 

Thus we conclude that mildly eccentric binaries could have been generated 
by fly-by encounters with low eccentricity pulsars below the line of "timing sen- 
sitivity limit" (group III pulsars). Some of the shorter P^rb binaries would again 
be circularised by gravitational radiation (see tgr contours in Fig [1] calculated 
using the formalism of Peters &: Mathews 1963). Group III binaries lie in the 
region tjiy < 10^^ years but they have not get eccentricities. Gravitational wave 
radiation can explain the circularity of tighter binaries among them, but others 
are either very young or lie far away from the cluster centre where vio/n^ has 
higher values. Extremely high eccentric binaries are probably results of either 
exchange or merger interactions. See Bagchi &: Ray (2009) for details. 
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